Optimizing the ionization and energy absorption of laser-irradiated clusters 



M. Kundu and D. Bauer 
Max- Planck- Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 

(Dated: February 2, 2008) 

It is known that rare-gas or metal clusters absorb incident laser energy very efficiently. However, 
due to the intricate dependencies on all the laser and cluster parameters it is difficult to predict 
under which circumstances ionization and energy absorption is optimal. With the help of three- 
dimensional particle-in-cell simulations of xenon clusters (up to 17256 atoms) we find that for a given 
laser pulse energy and cluster an optimum wavelength exists which corresponds to the approximate 
wavelength of the transient, linear Mie-resonance of the ionizing cluster at an early stage of negligible 
expansion. In a single ultrashort laser pulse, the linear resonance at this optimum wavelength yields 
much higher absorption efficiency than in the conventional, dual-pulse pump-probe set-up of linear 
resonance during cluster expansion. 

PACS numbers: 36.40.Gk, 52.25.Os, 52.50.Jm 



I. INTRODUCTION 

The interaction of rare-gas and metal clusters with in- 
tense laser light has drawn close attention during the last 
ten years. Reasons for this large interest are (among oth- 
ers) the high charge states [l|, [2, Jj, 11, M Js H, [a, 01 and the 
high energies of both ions [ u S M IiL 181,19, flfj, ITHH^ 
and electrons [H, E2> EE [la] observed (see [lrj] for a 
recent review). 

The laser-cluster interaction scenario may be qualita- 
tively summarized as follows. Electrons leave their "par- 
ent" atoms ("inner ionization") and absorb further laser 
energy while moving in the cluster potential formed by 
the ionic background. Some electrons may leave the clus- 
ter ("outer ionization"), leaving behind a net positively 
charged nanoplasma. The total electric field consists of 
the laser plus the space charge field. It may exceed the 
pure laser field in certain spatial regions and thus may 
liberate further electrons from their parent ions (which 
would remain bound if there was the laser field alone). 
This enhanced inner ionization is called "ionization igni- 
tion" [H Oil [H H3, HH Ionization ignition locally 
continues until the total field drops below the thresh- 
old field required to liberate further electrons. The ionic 
background expands because of Coulomb repulsion and 
hydrodynamic pressure, which ultimately leads to the en- 
ergetic ions observed in experiments. Because of the com- 
plex interplay between inner ionization, outer ionization, 
and ionic expansion it is far from trivial to predict quan- 
titatively mean or highest charge states, the absorbed 
laser energy, or other observables and their dependence 
on the laser parameters and the cluster decomposition. 

One of the goals in laser-cluster experiments and simu- 
lations is to convert as much laser energy as possible into 
energetic particles. This can be achieved by optimizing 
the outer ionization degree, i.e., by removing as many 
electrons as possible from the cluster in order to gener- 
ate high charge states so that the asymptotic ion energy 
(and thus the total absorbed energy) after Coulomb ex- 
plosion is largest. 

One way to increase the charge states and the ion en- 



ergy is to dope a cluster with atomic/molecular species 
of low ionization potential [H, HH . An almost two-fold 
increase of the highest charge states were obtained exper- 
imentally with argon clusters doped with water molecules 

Experimental results for xenon and silver clusters em- 
bedded in helium droplets were reported in Ref. [26[ . The 
pulse duration and the sign of the chirp of a laser-pulse 
also affect ion charge states and ion energies 27]. En- 
hanced inner ionization of rare-gas and metal clusters ir- 
radiated by a sequence of dual laser pulses were observed 
pdl l26l |28| . |29| experimentally. In these kinds of exper- 
iments one should adjust the delay time between pump 
and probe pulse such that the cluster expands sufficiently 
to meet the linear resonance WMic = with the probe 
pulse, where lumic is the Mie-plasmon frequency and u\ 
is the laser frequency. Vlasov simulations |26| and semi- 
classical simulations [3(| of a small Xe4o cluster subject 
to such a pump-probe setup showed an enhancement of 
the ion charge states. An optimum control multi-pulse 
simulation has also been performed [3l| . 

In this work we investigate the effect of the laser wave- 
length by three-dimensional particle-in-cell (PIC) simu- 
lations. The goal is to find an optimum wavelength for a 
fixed laser intensity and a given cluster. At this optimum 
wavelength (which turns out to be in the ultraviolet (UV) 
regime for the Xe clusters under consideration) a single 
ultrashort laser pulse is shown to be much more efficient 
than the "conventional" dual-pulse pump-probe setup. 

Experimental signatures of enhanced x-ray yields and 
high charge states at short wavelengths [H, [H, [34[ indi- 
cate a clear impact of the laser wavelength on the laser- 
cluster interaction. Free electron laser (FEL) cluster ex- 
periments [35| at the DESY facility, Hamburg, and a re- 
cent x-ray laser-cluster experiment (36j down to wave- 
lengths < 100 nm also showed enhanced ionization. 

Contrary to our findings recent molecular dynamics 
simulations [371 ] concluded that (i) there is no influence 
of the laser wavelength on the charging of clusters in the 
regime 100 — 800 nm for a laser intensity « 10 16 Wcm -2 
and (ii) that linear resonance plays no role, thus threat- 
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ening the basis of the nanoplasma model [38|. Simi- 
lar conclusions were reported by the same authors in 
Refs. [Hill. 

We consider short laser pulses in this work. Most of 
the earlier works were reported for the lo ng-p ulse reg ime 
where linear resonance (LR) absorption [2(| [3(| l3ll. l38l. 
EL E2, EH, El, Hf| occurs during the expanding phase of a 
cluster when the Mie-plasma frequency drops sufficiently 
so that the laser frequency can be met. At the time of 
LR the space charge field inside the cluster is strongly 
enhanced, leading to efficient ionization ignition. 

The paper is organized as follows. In Sec. [TT] we 
briefly describe the simulation method and discuss the 
ionization of a cluster by a short laser pulse in Sec. Mil 
In Sec. IIVI pump-probe simulation results are presented 
while Sec. |V] is devoted to the laser wavelength depen- 
dence of the cluster dynamics. A possibility to achieve 
100% outer ionization is also discussed in Sec. IVl before 
we summarize the work in Sec. I VII Unless stated other- 
wise we use atomic units (a.u.). 



II. DETAILS OF THE SIMULATION 

Details of our PIC code are already described in 
Refs. 0,113. 

For the inner ionization we apply the so-called Bethe- 
rule or over-the-barrier ionization (OBI) model (48[. Ac- 
cording to OBI an atom or ion is ionized if the total field 
satisfies 



z (R 3 ,t)\>I*(Z)/4Z 



(1) 



at individual ion locations Rj. Here Z is the charge 
state (after ionization), I p (Z) is the respective ionization 
potential and E sc (Rj, t) is the space charge field. In the 
absence of E sc (Rj , t) ionization is caused by the laser 
field, known as optical field ionization (OFI). 

In this work we shall vary the wavelength down to 
100 nm at an intensity 5 x 10 16 Wcm -2 , which raises 
the questions (i) whether such lasers are available and 
(ii) whether the Bethe-rule ([T]) is applicable. With the 
development of new generation FEL lasers [49|, [5(| the 
answer to (i) is clearly affirmative. As regards the ion- 
ization model (ii), at short wavelengths ionization rather 
proceeds via multiphoton ionization than via tunneling 
or over-the-barrier ionization so that the Bethe-modcl 
(where the ionization probability switches from zero to 
unity once a certain threshold field is reached) may not 
yield the precise charging dynamics of the clusters at 
short wavelengths. However, the final charge state dis- 
tribution should remain unaffected by the details of the 
ionization model [40| at least qualitatively. 

Only collisionless absorption mechanisms are incorpo- 
rated in standard PIC simulations. The neglect of par- 
ticle collisions is an approximation which is the more 
valid the smaller the clusters, the higher the laser in- 
tensities, and the lon ger the laser wavelengths are (see, 
e.g., [H, S3, HI S3, 1m HI)- Since collisions can be 
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FIG. 1: (Color online) Maximum ion charge -Z max (thick 
solid), minimum ion charge 2 m in (thin solid), average ion 
charge ,Z av (bold black) and the ion charge predicted by 
OFI alone (dashed) vs peak laser intensity for (a) a Xe2i76 
cluster of radius Ro m 3.54 nm and (b) a Xei7256 clus- 
ter of radius Ro ~ 7 nm in an n — 8 cycle laser pulse 
E\(t) = Eq sin 2 (ui\t/2n) cos(u>\t) of wavelength Ai = 800 nm. 



expected to increase both energy absorption and ioniza- 
tion, our results from collisionless PIC calculations may 
be considered close to reality for wavelengths ^ 400 nm 
but slightly underestimating the real charge states and 
the real absorbed energy for shorter wavelengths. 



III. IONIZATION OF A CLUSTER BY A 
SINGLE SHORT PULSE 

First we study the response of a xenon cluster in a 
linearly polarized n = 8-cycle laser pulse of electric field 
strength E\(t) = Eq sin 2 (wit/2n) cos(wit) and wavelength 
Aj = 800 nm. Different ionic charge states are self- 
consistently produced during the laser pulses according 
to the Bethe rule tjlj. 

Figure []Ji shows the maximum charge state 2 max , the 
minimum charge state Z m i n , and the average charge state 
Z av (defined as the total charge of the cluster divided by 
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the number of atoms N) and the charge state predicted 
by the OFI ( "Z-Bethe" curve) vs peak laser intensity for 
a XeAr cluster (N = 2176) of initial radius i?o ~ 3.54 nm 
after the pulse (i.e., after 22 fs). The maximum charge 
state Z max varies from Z = 5 to Z — 26 as the laser inten- 
sity increases from 2.5 x 10 14 Wcm -2 to 7.5 x 10 17 Wcm~ 2 . 
The higher value of Z niaK above the value predicted by 
the OFI is clearly due to ionization ignition. Those max- 
imum charge states i? max are mainly acquired by the 
ions at the cluster periphery where the space charge field 
is highest. Inside the cluster the total field falls be- 
low the ionization thresholds due to the decreasing space 
charge produced by the ionic background as well as due 
to the screening of the laser field by the cluster elec- 
trons. The ions close to the cluster center have minimum 
charge states Z m { n = 2 — 20 at laser intensities between 
2.5 x 10 14 Wcm- 2 - 7.5 x 10 17 Wc m - 2 . The value of Z min 
remains much lower than predicted by the OFI for almost 
all laser intensities < 5.0 x 10 17 Wcm~ 2 . The average 
charge Z av remains close to (but slightly higher than) the 
OFI predicted values at intensities < 7.5 x 10 15 Wcm~ 2 . 
Also Z max = 8 and Z m - m = 3 do not change between 
the intensities 10 15 Wcm -2 - 7.5 x 10 15 Wcm" 2 but Z av 
increases slowly as more ions from the cluster center to- 
wards the periphery acquire higher charge states 3^8. 
The value of 2 max remains constant, Z = 8, due to 
the removal of all electrons from the 5s 2 p 6 shell of the 
Xe atoms close to the cluster boundary. As the inten- 
sity s» 7.5 x 10 15 Wcm~ 2 is approached the laser field is 
strongly shielded from the central part of the cluster, and 
outer ionization as well as ionization ignition tend to sat- 
urate. As a consequence Z av grows slowly between the 
intensities « 5 x 10 15 Wcm~ 2 — 10 16 Wcm~ 2 . Unless a 
threshold intensity « 10 16 Wcm~ 2 is crossed further elec- 
trons from the cluster cannot be removed, which was al- 
ready seen in previous model and numerical calculations 
[IE [Hi [HI [55(]. At higher intensities > 10 16 Wcm~ 2 outer 
ionization and ionization ignition increases again, leading 
to an increase of Z av beyond the values predicted by the 
OFI due to the strong increase of both Z max and Z m i n . 

It is commonly believed that ionization ignition be- 
comes increasingly pronounced with increasing cluster 
size. Figure [TJd shows Z max , Z m ; n , Z av , and the charge 
states predicted by the OFI vs the peak laser intensities 
for a bigger Xejv cluster (N — 17256) of initial radius 
i?o ~ 7 nm. Z max varies between 8 — 26, exceeding again 
the charge states predicted by OFI alone. Below the in- 
tensity 10 17 Wcm~ 2 Z max is higher by a factor of w 2 
compared to the OFI value ("Z-Bethe" curve). Although 
2 max remains much higher, the average ion charge Z av (in 
Fig. Hfc) is below the charge states predicted by the OFI 
for most of the laser intensities. Most of the ions closer 
to the cluster center acquire charge states Z m ; n = 2 — 10 
which are even lower than for the smaller cluster (Fig.[T|a,) 
at the corresponding intensities. Hence, ionization igni- 
tion is indeed responsible for the highest charge states 
Z max which increase with the cluster size (as seen in 
Fig. []}. However, exactly because of the same mecha- 



nism a bigger cluster will capture more electrons (whose 
outer ionization would require much higher laser intensi- 
ties than in the case of a smaller cluster) . The presence of 
more electrons in the central region will screen the laser 
field more efficiently. As a result both Z av as well as Z m i n 
(in Fig. [TJd) drop below the corresponding values for the 
smaller cluster (Fig. [TJi) . 

We conclude that an increasing cluster size (and thus 
increased ionization ignition of, at least, the ions located 
close to the cluster boundary) does not always lead to 
a higher average charge state. Our aim is to increase 
not only the highest charge states but also the average 
ion charge beyond the OFI predicted value through the 
charging of more ions in the central part of the cluster. 
In the following sections we study several approaches to 
achieve this goal. 



IV. IONIZATION BY DELAYED PULSES: A 
PUMP-PROBE SIMULATION 

In this section we illustrate the "pump-probe" method 
frequently employed in laser-cluster experiments. In this 
method an initial pump-pulse ionizes the cluster. The 
cluster expands freely before, after a delay time, a probe- 
pulse hits the expanding cluster. The interaction of this 
probe pulse with the cluster will sensitively depend on 
the cluster size and thus on the delay time. We revisit 
such a scenario in our current work since it will allow us 
to compare the efficiency of laser energy absorption for 
such a standard pump-probe method with the single UV 
pulse scenario which will be introduced in Sec. fVl 

The laser field profile is of the form E\(t) = 
E sin 2 (7rf/nT) cos(wit) for both pump and probe pulse. 
The time period T is chosen with respect to the wave- 
length 800 nm. For, say, n — 4 the product nT deter- 
mines the total pulse duration ss 11 fs. The pulse enve- 
lope and intensity are kept the same for all cases under 
study, i.e., the laser energy in all pulses is the same too. 

Figure [2Ji shows the results for the Xei7256 cluster of 
initial radius Rq w 7 nm at an intensity 5 x 10 16 Wcm -2 
when both the pump and the probe pulse have the same 
wavelength 400 nm. The average charge Z av , the scaled 
Mic-frequency WMie/wi, the total absorbed energy .Etot 
(electrostatic field energy plus the kinetic energy of elec- 
trons and ions), the normalized cluster radius R{£)/Rq, 
and the laser fields are plotted vs time (in units of the 
period T) . During the first four laser cycles of the pump- 
pulse the average charge state rises to Z av w 11, the 
frequency WMie/^i ~ 2.5 and E to t ~ 2.0 x 10 7 while 
the cluster expansion is insignificant. The total energy 
E t ot ~ 2.0 x 10 7 corresponds to the average energy ab- 
sorbed per ion E to t/N 31.4 keV. After the pump-pulse 
the cluster evolves freely and Z av , E to t remain unchanged 
but c^Mie/wi drops due to the expansion. Note that the 
cluster radius R(t) (defining the cluster boundary) cor- 
responds to the distance of the most energetic ions from 
the cluster center. At the boundary, however, the cluster 
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FIG. 2: (Color online) Average ion charge 2 av , scaled Mie- 
frequency ujmic/uJi, laser field E\ (in atomic units), normal- 
ized cluster expansion radius R/Ro and total absorbed en- 
ergy _Etot (in atomic units) vs time (in 800 nm laser cycles) 
for a Xei7256 cluster of radius Rq ~ 7 nm. The peak in- 
tensity 5 x 10 16 Wcm -2 is the same for (a) pump of wave- 
length 400 nm (probe, 400 nm) and (b) pump of wavelength 
400 nm (probe, 200 nm). The laser field is of the form 
E\(t) — Eo sin 2 (irt/nT) cos(wii) with n = 4 and one laser 
cycle T corresponding to the wavelength 800 nm. A minor 
increase in iJtot after the pulses is an artifact of PIC simula- 
tions. 



potential is anharmonic. Hence using R(t) for the cal- 



culation of the Mie-frequency UMic{t) — y/N£ av /R 3 (t) 
the latter is und erestimate d. Instead we use the defini- 
tion ujuie(t) = \/Qb(t)/Ro (as in Ref. [13]) where Q\,{t) 
is the total ionic charge within the initial cluster radius 
Ro where the cluster potential is close to harmonic at all 
times. 

After 44 laser cycles WMic approaches the linear res- 
onance (dashed line) with respect to the fundamental 
400 nm, i.e., WMie/^i = 1- The probe pulse of wave- 
length 400 nm is applied with a delay of i=s 42 laser cycles 
such that the peak of the pulse approximately coincides 
with the resonance time. Due to the linear resonance the 
average charge and the absorbed energy rises abruptly 



up to the value Z av = 14 and E to t ~ 7.5 x 10 7 , respec- 
tively. Such a pump-probe simulation clearly illustrates 
that the linear resonance indeed plays a role in the clus- 
ter dynamics. More energy is absorbed, leading to higher 
charge states. These results are i n agree ment with hydro- 

However, linear 



dynamic and Vlasov simulations [111 
resonance is met only after a relatively long time when 
the cluster has already expanded significantly (as seen 
R(t)/Ro w 8 in Fig. [2^,). Ionization ignition and laser 
energy absorption in such a low density plasma is ex- 
pected to be less efficient compared to the case where 
linear resonance occurs before the cluster expands signif- 
icantly. 

While keeping the 400 nm pump as above we now as- 
sume a probe wavelength of 200 nm for the purpose of hit- 
ting the linear resonance at an earlier time when the clus- 
ter is more compact. The energy in the probe pulse is the 
same as in Fig. [2ji. Figure [2]d shows the result analogous 
to Fig. [2^,. The average charge and the absorbed energy 
now increase up to Z av = 18.5 and E tot ~ 9 x 10 7 which 
are higher than in Fig.[2Ji after the probe. With the pulse 
energies being the same in both cases a higher efficiency 
of energy absorption in the second scheme (Fig. is 
obvious. The reason is the smaller cluster size at the 
time of linear resonance (R(t)/Ro < 1.5) and the higher 
space charge field related to it. Similar findings from ex- 
periments have been reported in Ref. [Ill ] . In passing we 
note that the average charge Z av w 11 in Fig. [2] due to 
the pump (at 400 nm) exceeds Z av « 8 in Fig. [TJd (at 
800 nm) for the same cluster and the same laser inten- 
sity 5.0 x 10 16 Wcm~ 2 despite the higher pulse-energy in 
Fig. [TJb because of the twice longer pulse. 

In the following section we study the wavelength de- 
pendence of the average charge states and the laser en- 
ergy absorption. 



IONIZATION AT DIFFERENT 
WAVELENGTHS 



Does the average charge state and the absorbed en- 
ergy for a given cluster increases with decreasing laser 
wavelength? One could expect that for a certain wave- 
length the linear resonance during the initial ionization 
stage when the Mie-frequency rises from zero to its max- 
imum value becomes important. For long wavelengths 
this early resonance is passed so quickly due to the rapid 
charging of the cluster that any indication of a resonance 
is washed out. 

We assume the same laser field profile E\(t) = 
E sin 2 (Trt/nT) cos(wit) as in SeaHTDwith the same pulse 
duration, pulse energy, and intensity 5.0 x 10 16 Wcm -2 
so that the number of laser cycles in the pulse depends 
on the wavelength. The laser wavelength is varied in the 
range 800 — 100 nm. Note that in the following we spec- 
ify times and pulse durations in units of laser periods at 
800 nm (corresponding to T « 2.6 fs). 

Figure [3] shows the average charge state Z av and to- 
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FIG. 3: (Color online) Average ion charge (a) and total 
absorbed energy Etot (b) vs laser wavelength after 1, 2, 3 and 
4 laser cycles (at 800 nm) for a Xei.7256 cluster of radius -Ro « 
7 nm. Other parameters as in Fig. [2] 



FIG. 4: (Color online) Normalized cluster expansion radius 
R/Ro (a) and outer ionization degree r\ (b) vs laser wave- 
length, corresponding to Fig. [3] 



tal absorbed energy E toi vs the laser wavelength for the 
Xei7256 cluster of radius i? ~ 7 nm after t — 1,2,3,4- 
laser cycles at 800 nm. The value of Z av increases in 
time (in Fig. [3^,) for all wavelengths. Ionization mostly 
occurs before t = 2 cycles when the peak of the pulse is 
reached. After that the space charge field is high enough 
to generate further charge states between 2 — 3 cycles. 
Z av does not change anymore between 3 — 4 cycles, in- 
dicating a saturation of inner ionization. The average 
charge state Z av increases from Z = 8 to a maximum 
value Z aw ss 25 as the laser wavelength is decreased from 
the infrared 800 nm down to the UV wavelength 125 nm. 
It means that the sub-shells 4s 2 p 6 d 10 5s 2 p 6 of almost all 
atoms are empty at 125 nm. A further decrease of the 
wavelength causes Z av to decrease gradually to a smaller 
value Z av w 9 at 50 nm. 

Figure shows a similar qualitative behavior of the 
absorbed energy both in the time domain and in the 
wavelength domain. The energy E to t is maximum at 
the same wavelength Ai = 125 nm. Although the laser- 
pulse energy is the same in all cases the increased ab- 



sorption at 125 nm, leading to a marked increase of the 
average charge up to a value Z av w 25 clearly shows 
that wavelength effects are undoubtedly important. One 
may compare the absorbed energy and the average charge 
with the dual-pulse simulation results in Fig. [2j The ab- 
sorbed energy E to t ~ 34 x 10 7 a.u. and the average charge 
Z fa 25 are much higher in the present case around the 
laser wavelength 125 nm compared to the respective val- 
ues E tot « 9 x 10 7 and Z av « 18.5 in Fig. \2k>- The 
absorption is w 3.78 times higher than in Fig. [^b- More- 
over, in the dual-pulse case the total laser-pulse energy 
was twice higher. Therefore, the absorption efficiency is 
augmented further by a factor of two. 

In Fig. [¥] we plot the normalized expansion radius 
R(t)/Ro (Fig. |4ji) and the outer ionization degree ij 
(number of total electrons outside R(t) divided by the 
total number of electrons produced, NZ av (t), in Fig-dJa) 
vs the laser wavelength corresponding to the results in 
Fig. [3l The radius R(t) and the outer ionization de- 
gree 77(f) go hand in hand with the absorbed energy 
Etot{t) and the charge Z av (t). After four cycles the 
cluster has expanded very little, R{AT)/Rq w 1.225 at 
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FIG. 5: (Color online) Scaled Mie-frequency uiMie/ui vs time 
(in 800 nm cycles) for wavelengths Ai = 800 — 112.5 nm and 
the laser and cluster parameters of Fig. [3] 



Ai « 125 nm, although the average charge Z av w 25 is 
very high compared to Fig. O With such an insignifi- 
cant expansion the space charge field can be considered 
optimized, leading to maximum ionization ignition. The 
ignition field (i.e., the space charge field due to the ionic 
background) under the assumption that all electrons are 
removed reads Ei g (t) « N Z av (t) / 'R(t) 2 . Using R(t) from 
Fig. QJi and Z av (t) from Fig. [5^, one obtains at 125 nm 
E ig (2T) « 20.0, E ig {3T) « 21.Q and E ig (4T) w 16.2, if all 
electrons are removed (i.e., 77 = 1). The expected ignition 
field 21 a.u. is maximum near the pulse peak around 
2 — 3-cycles, thereafter decreases to E ig (4T) w 16.2 due 
to an expansion R(4T)/Rq w 1.225 and no further cre- 
ation of charge states. Note that the peak laser field is 
only Eq rs 1.19. Therefore the enhanced ionization is cer- 
tainly due to the ignition field. However, at Ai rs 125 nm 
77 w 0.6 in Fig.QJ), meaning that 40% of the electrons are 
still inside the cluster. The presence of these electrons 
lowers Ei g compared to the above ideal case of i] — 100% 
outer ionization, and one may argue that Ei g is not yet 
optimized. However, even if r\ = 100% outer ionization is 
achieved for the above laser field intensity the maximum 
total field is w 22 which is still insufficient to produce 
a higher average charge Z av = 27 (requiring a threshold 
field gl 24 according to OFI). Hence the average charge 
state is optimized. This will be shown explicitly at the 
end of this section where we actually achieve r\ w 100% 
for this cluster. 

The above results clearly show that there exists a cer- 
tain wavelength at which the laser-cluster coupling is very 
efficient. Such a nonlinear dependence of the absorbed 
energy and average charge state on the laser wavelength 
indicates a resonance around 125 nm in Fig. [3] and |4] 
To investigate this further, we plot in Fig. [5] the scaled 
Mie-frequency LOy\\ a {€) j lo\ vs time. The dashed line indi- 



cates the linear resonance. Charging of the cluster starts 
around 0.3 cycles for all wavelengths by OFI, leading 
to an abrupt increase of WMie(*)/wi for the longer wave- 
lengths while for the shorter ones the increase proceeds 
slower. As a result the plasma is overdense during the en- 
tire pulse for the long wavelengths but stays close to the 
linear resonance for the shorter wavelengths. The more 
time is spent near the linear resonance, the higher is the 
energy absorption and the average charge state, as seen 
in Fig. [3] At the wavelength 125 nm the resonance is met 
at the peak of the pulse so that the energy absorption is 
particularly efficient. 

We now discuss the time evolution of the space charge 
field E® c (along the laser polarization) at different posi- 
tions inside the cluster to further illustrate the resonance 
at short wavelengths, leading to efficient ionization igni- 
tion and the generation of high charge states. Figure [5] 
shows the space charge field E* c and the laser field E\ 
vs time at radial positions 0.24i?o, 0.48i?o, 0.72i?o and 
0.96i?o for four different wavelengths. Figure [7] shows the 
corresponding phases of Ef c with respect to the driving 
laser field. 

At the long wavelength 800 nm, E* c inside the cluster 
at radii 0.24i? , 0.48i? , 0.72i? mostly oscillates with 
a phase 5 rs tt while S ~ at the boundary (i.e., at 
0.96i?o)- This is clearly what one expects from an over- 
dense plasma: screening of the laser field in the cluster 
interior but an opposite behavior outside the electron 
cloud. The oscillation of the space charge field arises due 
to the oscillations of the electrons trapped inside the clus- 
ter. These electrons form approximately a sphere which 
is smaller than the cluster due to outer ionization. If 
the electron cloud was rigid and did not cross the clus- 
ter boundary the phase should be exactly tt and inside 
and outside, respectively, if the plasma is overdense, and 
opposite in the underdense case. In reality, the bound 
electron population changes and the electron sphere is 
neither rigid nor has it a sharp boundary, resulting in 
phase distortions and deviations from the idealized case, 
as seen in Fig. [7^,. 

Figure [5] confirms explicitly that the total field at the 
boundary is highest and therefore leads to the highest 
ionic charge states while E* c {t) almost nullifies the laser 
field in the strongly overdense regime. The maximum 
value of the total field at the peak of the pulse is rj 4.0 a.u. 
(x-component only) which is sufficient to produce charge 
states up to Z w 18 (also seen in Fig. QJ>). An additional 
contribution (up to a factor y3) to the total field comes 
from the y and ^-components of the space charge field. 

At 200 nm the amplitude of E* c around t = 1.5 cycles 
at 0.24i?o increases up to 5 a.u. which, after addition 
to the laser field, is sufficient to produce charge states 
Z = 18 even inside the cluster. After w 1.75 cycles E* c 
at 0.72i?o behaves similarly to that at 0.96i?o, i-e., the 
laser and space charge fields at 0.96i?o and 0.72i?o are 
now approximately in phase. Figure [SJs shows that the 
total field w 4 — 8 between 0.96i?o and Q.72R , producing 
charge states Z w 18 — 23. However, due to the screen- 
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FIG. 6: (Color online) The ^-component of the space charge 
field at radial distances 0.247? , 0.487? , 0.72i? , and 0.96i? 
inside a Xei7256 cluster of radius Ro w 7 nm and the laser 
field E\(t) = Eo sin 2 (ivt/nT) cos(oJii) of peak intensity 5 x 
f0 16 Wcm~ 2 vs time (in periods corresponding to 800 nm) at 



(a) Ai = 800, (b) 
112.5 nm. 



Xi = 200, (c) Ai = 125, and (d) Ai = 



FIG. 7: (Color online) Phase of the space charge field with 
respect to the laser field at different radial distances corre- 
sponding to Fig. [6]vs time corresponding to Fig. [6] 
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ing of the laser field inside the cluster, many atoms there 
have only charge states Z < 18 so that the average charge 
state is Z av w 18 in Fig. [3] From Fig. [7)d it is seen that 
before t — 1.75 cycles the phase of E* c at 0.72i?o ap- 
proximately follows the phase at the smaller radii 0.24i?o 
and 0.48i?o since the plasma is evolving from under to 
overdense. Then, with increasing outer ionization and 
thus shrinking electron sphere, E* c at 0.72i?o drops and 
approaches the behavior for 0.96i?o- 

At the resonant wavelength 125 nm violent oscillations 
of the electron cloud are driven, leading to a particularly 
high total field everywhere inside the cluster and an aver- 
age charge state Z av s» 25 in Fig. [3] Higher charge states 
Z > 26 are not produced because of the high thresh- 
old field w 24 necessary to crack the M-shell. One may 
argue that the presence of 40% electrons inside the clus- 
ter (in Fig. 0b) will deplete the field inside significantly. 
However, one should keep in mind that at resonance the 
electron cloud oscillates with a large excursion, exposing 
a substantial part of the naked ionic background, lead- 
ing to an enhanced "dynamical ionization ignition" [l8j 
which can produce higher charge states than expected 
from the laser field alone even inside the cluster. Finally, 
after t = 3 cycles E* c at the boundary drops due to the 
cluster expansion. As expected, the phases plotted in 
Fig- Eh fluctuate around tt/2 throughout the cluster once 
the resonance condition is met. 

At 112.5 nm the plasma remains underdense. Fig- 
ure [5Ji shows that the space charge field amplitudes drop 
compared to those in Fig. [St, yielding less ionization ig- 
nition and absorbed energy, similar to the 200nm-case. 

Even in the optimal 125 nm case presented so far only 
60% of the generated electrons were removed from the 
cluster (visible in Fig. UJd). Therefore outer ionization 
and ionization ignition was certainly not optimized. We 
argued that even if the remaining 40% electrons were 
removed, the average charge state would not be signifi- 
cantly increased as compared to that shown in Fig.[3^. To 
prove that, we performed PIC simulations for the same 
cluster and the same peak intensity but now employing 
two consecutive pulses (shown in Fig. [5J . The first pulse 
of resonant wavelength 125 nm with respect to the still 
compact cluster is ramped up over four 800 nm-cycles 
and held constant afterwards up to 8 cycles (the details 
of how the pulse is ramped down do not matter; therefore 
it is simply switched off abruptly) . At t = 8 cycles a sec- 
ond pulse is switched on (over 2 cycles) whose frequency 
is resonant with the Mie-frequency around t = 10 cycles. 

After the first pulse the cluster doubled its radius, and 
the outer ionization degree amounts to r/(R) w 0.8 so 
that 20% electrons are still inside the cluster of radius 
R(t) while w 10% are inside a sphere of radius Ro, The 
average charge Z av in Fig. [8p does not change signifi- 
cantly compared to Fig. although the pulse energy 

per unit area J^ T £f(i) dt is » 3 — 4 times higher. 

The purpose of the second pulse shown in Fig. [5b is 
the removal of the residual electrons. Although almost 
95% outer ionization within the expanding radius R and 
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FIG. 8: (Color online) Variation of (a) normalized cluster 
radius R/Ro, degree of outer ionization r](R),r)(Ro) within 
R and Ro, laser field E\ (in atomic units) (b) average ion 
charge iJ a v and total absorbed energy £tot in time (in units 
of the laser period T corresponding to 800 nm) for the Xei7256 
cluster of Fig. O The two laser pulses of wavelength 125 and 
415 nm, respectively, are included in (a). The peak intensity 
is the same as in Fig. [3] for both pulses. 



99% within Rq are achieved, no higher charge states are 
created. The absorbed energy also does not rise signifi- 
cantly so that the higher input energy invested into the 
two pulses does not pay off. Hence a single, short UV- 
pulse of wavelength 125 nm turns out to be optimal with 
respect to fractional energy absorption and generation of 
a high average charge state under the conditions consid- 
ered. 



VI. SUMMARY 

In summary, we studied the interaction of xenon 
clusters with intense short laser pulses using a three- 
dimensional PIC code. Our aim was to optimize for a 
given cluster the laser energy absorption and the gen- 
eration of high average charge states. The latter will 
then lead to energetic ions upon Coulomb explosion. We 
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showed that for a given laser intensity an optimal laser 
wavelength exists that, under the typical conditions stud- 
ied in this work, lies in the UV regime. Energy absorp- 
tion is optimized when resonance is met during an early 
stage of the dynamics when the cluster is still compact. 
The conventional, long-pulse linear resonance during the 
expansion of the cluster is less efficient. 
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